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Recent experimental observations on the 'bulk' features of particle production at high 
(pseudo)rapidities will be reviewed. This kinematic region is of interest mostly because 
of its relevance to the theoretical description of initial state effects of nuclei at ultra- 
relativistic energies. Measurements of the charged hadron multiplicity density as well 
as the pseudorapidity dependence of the elliptic and directed flow exhibit a remarkable 
scaling property as a function of collision energy. This scaling seems to hold for pions and 
even photons and J/^-s, but is violated for protons. The special role of baryons will be 
discussed using selected results on nuclear transparency and baryon stopping. 

1. INTRODUCTION 
1.1. Motivation 

The focus of research in relativistic heavy ion collisions initially has been toward the 
mid-rapidity region. This is justified in the context of the search for the new phase of 
QCD at high temperatures (energy density), the Quark Gluon Plasma, the emergence of 
which was expected at midrapidity. In addition, an extended Lorentz-invariant central 
plateau was predicted in the (charged) particle rapidity distribution [ [T] . Most experi- 
ments have been designed with good midrapidity-coverage in mind. At earlier fixed target 
experiments, particle detection closer to beam (target) rapidity was less difficult than at 
high energy colliders. 

In heavy ion collisions, the existence of a boost-invariant central region is not confirmed 
experimentally. On the other hand, several observables in the high rapidity region are 
found to scale in a simple way as a function of energy. Furthermore, most of the baryons 
are emitted at high rapidity, and their energy loss is of primary interest, since it is related 
to the created energy density. The total number of produced charged particles per par- 
ticipant nucleon (N part ) pair does not depend on the centrality of the heavy ion collision. 
There is a remarkable balance of the centrality dependent midrapidity and high rapidity 
particle yields [|2j, suggesting a deeper connection between these two regions. 

At RHIC, a suppression of high-p T particles at high rapidities is observed in d+Au colli- 
sions [E], which is interpreted as the saturation of the gluon density of the ultra-relativistic 
Au nucleus at low Bjorken-x values [Ej. Measurements at high rapidity are necessary to 
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Figure 1. Constant y (dashed lines, for proton mass), r\ (thin solid lines) and XBj (dash- 
dotted lines) contours in the pr vs. pz plane for the highest RHIC energy. (See the text 
for details.) The PHOBOS acceptance for very low pt [E] (low-XBj and low Q 2 ), and 
part of the BRAHMS acceptance at high rj (low-a^j and higher Q 2 ) are shown by ellipses. 



reach small enough x values but still have a large enough momentum transfer (Q 2 ) for this 
interpretation to apply. However, it is not easy to disentangle which physical processes 
(Cronin-effect, shadowing, energy loss, initial state gluon saturation, fragmentation of 
partons at high rapidities, etc.) influence the observed final hadron yields. 

There is a large amount of data available on collisions of simpler systems (p+p, P+A) 
partially still under analysis. Baryon stopping and nuclear transparency is a large field of 
study in itself, and is also related to basic questions about the carriers of baryon number. 
Parton energy loss in the QGP and in cold nuclei are at the center of attention ever 
since the discovery of the suppression of high pt particle yields ('jet quenching') [E]. 
The observed baryon 'anomaly' (unexpectedly high p/n + and p/ii~ ratios at high pt) [IE] 
suggests another connection between dynamics of baryons at midrapidity in the transverse 
direction, and baryon transport in the longitudinal direction A non-exhaustive set of 
measurements will be discussed here, to illustrate some of the thoughts outlined above. 

1.2. Kinematics of saturation physics 

Various kinematic variables are used to study high rapidity phenomena and saturation 
physics, where the largest depletion of gluon density is expected at low xej and high Q 2 ; 
experiments have to focus on that region. Fig. ^ shows constant rapidity (y, dashed lines, 
for proton mass), and pseudorapidity (rj, thin solid lines) lines of final state particles in 
the transverse (pt) vs. longitudinal momentum (pz) plane for the highest RHIC energy 
(100 GeV/nucleon colliding beams). As an illustration, constant Bjorken-x (xbj) lines of 
partons probed in the Au nucleus are also plotted (dash-dotted lines, ignoring mass) for 
a d+Au collision (with the Au nucleus moving from right to left) 1 . A fourth variable, 
Feynman-x (not shown) is defined as = 2pz / ^s NN , thus ip ~ XBj for large XBy 



1 We use a great simplification here: we think of the d projectile as the incoming 'lepton' and the final 
state particle as the outgoing 'lepton', and borrow the definition of xb} from Deep Inelastic Scattering. 
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Figure 2. Pseudorapidity density distributions of charged particles emitted in p(p) +p (left 
panel) and d+Au, p+Emulsion (right panel) collisions at various energies as a function of 
the rf — r] — Ubeam variable [|T7l [THJ [THJ I2DI- For details of the normalization of the data 
in the right panel see [I19j. 



1.3. Experimental capabilities 

At RHIC, several detectors cover the high rapidity region. The STAR experiment has 
the Photon Multiplicity Detector [|S] (2.3 < rj < 3.8), Forward TPC [Ej (2.5 < \rj\ < 4.0), 
Forward 7i~° Detector (3.3 < rj < 4.1) and the planned Forward Meson Spectrometer 
(2.5 < 7] < 4.0). The PHENIX experiment can use the muon stations to measure hadron 
yields at high rapidities (1.4 < |^| < 2.2), using the decay length dependence of the 
K ± — > ij ± i' i1 {V^) and 7r ± — > u ± i'^{V^) decays, or observing punch-through hardons [ 
ITT] . The BRAHMS spectrometers cover a large pseudorapidity region (77 < 3.5) with 
particle-ID using Ring Imaging Cerenkov Detectors [U2]. Finally, the largest 77 acceptance 
(although with no particle identification) is offered by the PHOBOS Multiplicity Array 
\r]\ < 5.4 [HH]. As examples, the PHOBOS acceptance for very low pt [Ej (low-^Bj and 
low Q 2 ), and the BRAHMS acceptance at high rj (Iow-^bj and higher Q 2 ) are illustrated 
in Fig. ^ The low xej and high Q 2 region is at high rj and high px- At RHIC, this 
kinematic window is rather narrow, however, a large momentum space window to study 
low-^Bj physics will be available, even at midrapidity, to the future LHC experiments. 
The CMS experiment plans to cover an extensive range of pseudorapidity (|?7| < 6.9) with 
its calorimeters [ITB]. 



2. ENERGY DEPENDENCE OF PARTICLE PRODUCTION AT HIGH r] 

2.1. Inclusive charged particle spectra 

One of the simple ways to study particle production at high rapidity is to compare 
charged particle r\ distributions measured at various collision energies. The r]'=r]-y^ eam 
variable is used. It is measurable without particle identification and approximately trans- 
forms the distributions to the rest frame of one of the colliding particles/ions. 

In proton-proton collisions, the so called limiting fragmentation [ 16J is seen from a few 
GeV to 900 GeV collision energy (e.g. at CDF, UA5 and ISR [|T7j). The left panel of Fig. 



4 



G. I. Veres 



CM 

10 
Q. 

■o 

1 



- I ' M I ' M 

1 ^200 GeV : 




ni30GeV 




; #62.4 GeV 




^19.6 GeV 




a) Au+Au 0-6% %^ 


b) Au+Au 35-40%^%^ 


PHOBOS , i^asfc 


i l i i t l i i i l t i i I i i i I ' 



s. o 



10 



a 200 GeV 
□ 130 GeV 

• 62.4 GeV 

* 19.6 GeV 



'-6 -4 




-2 2 -6 -4 -2 2 -2 -1 

T|' T|' if 

Figure 3. Left panel: distributions of pseudorapidity density of charged particles emitted 
in a) central, b) peripheral Au+Au collisions at various energies as a function of rj' = 
V — Vbeam, normalized by the average number of participant nucleon pairs. Right panel: 
ratio of rj' distributions in peripheral and central Au+Au collisions, normalized to the 
number of participant pairs [I21j. 



|2] shows that the charged particle pseudorapidity density distributions are approximately 
independent of collision energy over a large range of rj', which grows with y'i. 

As the right panel of Fig. El shows, similar scaling was observed for p+A collisions 
between ^/s NN = 6.7 and 38.7 GeV [HBj, and more recently in d+Au collisions at y/s NN = 
200 GeV [ UHl E0j ■ The observed energy- independence of these distributions holds both 
for the projectile and for the target reference frames, i.e. for rj ± ybeam [ITS]- 

In heavy ion collisions at RHIC energies, the BRAHMS and PHOBOS experiments 
have collected a large and systematic dataset on dN/drj distributions, at y/s NN = 19.6, 
62.4, 130 and 200 GeV energies in Au+Au [ E2| , and at 62.4 and 200 GeV energies 
in Cu+Cu collisions [EI]- Instead of a sizeable plateau at midrapidity, similarly to p+p 
and p+A collisions, an extended longitudinal scaling region around rf — is observed, 
as shown in the left panel of Fig. El One can visualize this scaling by imagining that 
the energies of the two colliding beams are adjustable separately. If one holds the energy 
of the first beam constant while increasing the energy of the second beam (in which the 
low-x gluons become more and more saturated), the particle density measured at a fixed 
7] increases only up to a limit, given by the limiting curve in Fig. |21 where it saturates. 
This longitudinal scaling holds for Au+Au events with the same collision centrality (same 
number of participants), but the limiting curve is different for central and for peripheral 
events. The ratio of dN/drj' distributions in peripheral to central Au+Au collisions is 
found to be energy independent to a remarkable precision over an order of magnitude in 
^fs, as shown in the right panel of Fig. 01 

Longitudinal scaling is compatible with the Color Glass Condensate picture. The appli- 
cation of gluon saturation and standard parton fragmentation to the high rapidity region 
leads to such a limiting curve [ 23J. This prediction is shown in the left panel of Fig. HJ 
and is in good agreement with the data. 

Of course, the theoretical description breaks down at a certain rj' value, as the probed 
x-Qj in the projectile increases. The centrality dependence of the dN/drj distributions is 
also correctly described by the Color Glass Condensate model [l2*lj. 
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Figure 4. Left panel: dN/drf distributions of charged particles emitted in central Au+Au 
collisions at various energies as a function of rf = r\ — ?/bcam, normalized by the average 
number of participant nucleon pairs, together with CGC calculations [I2S]- Right panel: 
extrapolation to LHC energy based on linear limiting curve and logarithmic rise of midra- 
pidity density with yfs [ 26. 

The predictions of various model calculations for dN/dn at rj = in central Pb+Pb 
collisions at the LHC ( v /s^=5500 GeV) range between 1000 and 7000. Based on the 
simple assumptions that the limiting curve of the dN/drf distributions is linear (as seems 
to be the case at RHIC energies) and that the observed logarithmic increase of dN/drjl^o 
with y/s NN continues to be valid [ |2H] , one can estimate the central dN/ drj to be about 
1100 and the total number of charged particles created to be about 14000 at the LHC [ 
121)] , see the right panel of Fig. HJ 

2.2. Identified particles 

The scaling properties discussed so far are more surprising if one considers the different 
production mechanism of various particle species. The abundance of the lightest mesons 
indicate the amount of entropy created (in a thermodynamical picture). The amount of 
net baryons (B — B) however is a conserved quantity, and the amount of stopping of 
the initial baryons in the colliding nuclei gives rise to the energy necessary for meson 
production as well as baryon pair production. 

There are extensive measurements of the rapidity distributions of charged pions [ 123 
I2H EH]- protons and antiprotons [123 EH EH] i n heavy ion collisions. As shown in Fig. 
for collision energies between 2.63 and 200 GeV per nucleon pair, the charged pion 
production, observing it in the rest frame of one of the projectiles, shows a similar scaling 
behaviour as the charged hadrons as a whole. The recent measurement of the neutral pion 
yield at high rapidities at 62.4 GeV [E2] follows the scaling observed for the charged pions 
rather precisely. A similar scaling of net proton distributions is not expected because of 
baryon conservation and a widening of the rapidity gap with increasing energy. The right 
panel of Fig. EJis a compilation of dN/dy distributions of protons (below y/s NN =5 GeV) 
and net protons (at ^/s NN = 17.27 and 200 GeV). 




The natural question arising is whether this longitudinal scaling at high rapidity is a 
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Figure 5. Rapidity (y'—y-ybeam) distributions of 7r and n + particles and net protons 
(p — p) produced in central Au+Au (Pb+Pb) collisions at y/s NN = 2.63, 3.28, 3.84, 4.29 [ 
6.27, 7.62, 8.76, 12.32, 17.27 [EEIED! and 200 GeV [Ellfil, and vr°-s at 62.4 GeV [ 

E2|. 



fundamental feature of particle production, or a remarkable coincidence, considering that 
many effects (baryon stopping, gluon saturation, collective effects, parton fragmentation 
etc.) play a role at the same time. Naturally, the scaling behaviour cannot hold equally 
precisely for all charged particles and for only mesons, since the dynamics of baryon 
stopping and production is fundamentally different. This could perhaps explain why the 
limiting curve for charged particles appears to be centrality dependent, as opposed to 
photons - originating mostly from 7r° decays [13*2*]). 

2.3. Feynman-x scaling of various final states 

A frequently used way of quantifying longitudinal scaling and target size dependence 
of particle production in p+A collisions is to study the a exponent in the relation a p A = 
aoA a . This is done by measuring the <t p a production cross sections as a function of the 
target nuclear mass, A, and as a function of the Feynman-x variable. The left panel of 
Fig. |H] is a compilation of data showing various final states and collision energy in p+A 
collisions [EH]- This shows that the attenuation as a function of xp is remarkably uniform 
within the experimental uncertainties. A more recent measurement of J/ty mesons in 
d+Au collisions at ^/s NN =200 GeV [ ESI i s compared to lower energy data (right panel 
of Fig. IHJ) and also exhibits the xf scaling with energy. However, the same scaling does 
not hold if the same data is plotted as a function of the XBj variable (referring to the Au 
nucleus. This indicates that the way in which gluon saturation plays a role in the J/ty 
production is strongly energy dependent (middle panel of Fig. |UJ). 

In the following section more details about the centrality dependence of hadron pro- 
duction (and baryon stopping) at high rapidities are discussed. 

3. BARYON STOPPING AND NUCLEAR TRANSPARENCY 

It is apparent from the previous sections, that the dynamics of the energy loss by the 
initial baryons in a heavy ion collision strongly influences particle production (at high 
rapidities). Baryon energy loss in the hot medium created in heavy ion collisions (A+A) 
and in cold nuclear matter (p+A, d+A) has been studied and compared extensively (e.g. 
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Figure 6. Scaling of the production of $,7r + ,7r ,p, p, n, A, A, K®, S°, K + , K and J/\& 
particles as a function of xp (left and right panels) and x-q] in the target nucleus (middle 
panel) in p+A [El and d+Au [E3] collisions at v /s^7=19, 39 and 200 GeV. 



[EH ESI EH EZj)- In addition, baryons are observed to behave differently than mesons 
when studying the 'jet quenching' of different species, leading to the baryon 'anomaly' at 
RHIC [E]. 

Among the earlier puzzles on nuclear transparency is the independence of the A/ A 
ratio on the size of the target nucleus (left panel of Fig. EJ) seen in minimum bias p+A 
collisions at 300 GeV beam energy on fixed target [EH] • A complete attenuation in the cold 
nucleus, which would restrict particle production to the nuclear surface, may seem to be a 
plausible explanation. Preliminary data on A hyperons at RHIC [EHj show that this ratio 
is independent of centrality in d+Au collisions (second panel in Fig. |7J). However, the 
dN/dy yield of net lambdas (A— A) strongly depends on target size and centrality, showing 
clear evidence of more baryon stopping for larger targets and central collisions, shown in 
the two rightmost panels of Fig. [7| (200 GeV proton beam on fixed Au and S targets [ HU] 
and d+Au collisions at y / s^=200 GeV [[39J). Similar observations have been made for 
the p/p ratio and for the net protons (p — p), see for example [ I40j . Interestingly, in d+Au 
and Au+Au collisions at v /s^"=200 GeV energy, the p/p ratio does not, or only weakly 
depends on centrality at midrapidity [HI], while the amount of net protons is a strong (in 
fact linear) function of the number of participants in Au+Au collisions at a/s nn =62.4 and 
200 GeV [E21- Other studies of p+Pb and Pb+Pb collisions at SPS have also shown that 
while there is no complete attenuation for baryons, the net proton distributions have a 
very strong centrality dependence. Remarkably, the stopping is stronger in central p+Pb 
than central Pb+Pb collisions, even if samples with the same average number of binary 
nucleon-nucleon collisions per participant are compared [ ESI ESI EZ] • 

In central Au+Au collisions at ^/s~=200 GeV, the average rapidity loss of the incoming 
baryons was estimated to be about two units [ &\ , and a similar value was obtained 
earlier for minimum bias p+A collisions [H3j. Furthermore, baryon stopping (the rapidity 
distribution of p — p) was found not to depend on the size of the projectile nucleus as long 
as it was smaller than the target nucleus and the collision was central [ I40j . 

The examples above show that it is not only necessary to understand baryon stopping to 
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Figure 7. A/A ratio vs. xf from [EH] in p+A collisions (300 GeV protons on fixed target); 
the same ratio in d+Au collisions at v /s^=200 GeV from [[3H| (preliminary data); the 
net A rapidity distribution from [Ell; and in min. bias p+A collisions at the SPS [ I40j . 

separate it from other physically important phenomena. To reach the correct conclusions, 
one must consider the relation of high- and mid-rapidity physics, small and large colliding 
systems and the centrality dependence of the various observables, all at the same time. 
Studies of baryon stopping may also be able to address profound physical questions about 
the nature and carriers of baryon number. As an example, the HIJING BB 2.0 model 
including baryon junctions can describe the rapidity dependence of both the p/p ratio and 
the net baryon density in 200 GeV Au+Au collisions [Ej. 



4. AZIMUTHAL ASYMMETRY OF HIGH-?] PARTICLE PRODUCTION 

One of the important 'bulk features' of hadron production at high rapidities is the 
asymmetry of the azimuthal distribution of particles around the beam direction, relative 
to the reaction plane ('flow'). This asymmetry can be elliptical, characterized by the 
t>2 coefficient, as well as off-center or 'directed', characterized by the v± coefficient of its 
Fourier-expansion. This flow is a signature of collective, final state effects, as opposed to 
the gluon saturation at low xbj, which is truly a feature of the initial state in the collision. 
Thus, the saturation picture alone cannot account for the large v% values observed in A+A 
collisions at RHIC. Furthermore, there is no reason to expect that similar longitudinal 
scaling would exist for the flow observables, as it did for the dN / dvj distributions in that 
picture. In fact scaling is found to hold [ED El) as can be seen in Fig. |H1 both for the 
and the V\ variables, in Au+Au collisions from 19.6 to 200 GeV 2 . The elliptic flow data 
are well described by an analytical hydrodynamic model assuming ideal fluid dynamics 
[ I3H] ■ This agreement illustrates how important final state interactions are for particle 
distributions at high rapidities. 



5. SUMMARY 



In summary, an extensive amount of data on hadron production at high rapidities has 
been collected and is currently under analysis and interpretation. From these data only 



2 Note that the strong rj dependence of V2 was not observed by the BRAHMS experiment 
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Figure 8. The v-i measure of elliptic collective flow (left panel) and the v\ measure of 
directed flow (right panel) in Au+Au collisions at y/s = 19.6, 62, 130 and 200 GeV energy, 
as a function of rf = r\ — ybeam [ HH EE| • (The centrality is 0-40% of the total inelastic 
cross section.) Theoretical curves were made using the Buda-Lund model [ I46j . 



a few results could be selected for discussion. Longitudinal scaling of dN/drj and dN/dy 
distributions has been observed as a function of energy for pions and charged hadrons but 
not for protons. Scaling features, the centrality dependence and dN/dr] distributions are 
successfully described by initial state gluon saturation models. However, baryon stopping 
(the influence of the valence structure of the projectile) also contributes to the final state 
hadron yields at high rapidity, following complicated dynamics having energy, centrality 
and target size dependence. Remarkably, another important 'bulk feature' of particle 
production, the azimuthal asymmetry, which is not explainable by initial state effects 
alone, still shows the same kind of longitudinal scaling (energy independence) in heavy 
ion collisions. Theoretical progress in describing high rj (low-x) particle production in p+A 
and A+A collisions is rapidly progressing and generating a lot of interest. Nevertheless 
for a clean experimental isolation of initial state modifications of the incoming projectile 
one might have to turn to more differential measurements. Ideas to separate the various 
physical phenomena at high rapidities along with the new experiments at the LHC, will 
guide our theoretical understanding. 
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